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We study equilibration of quantum Hall edge states at integer filling factors, motivated by experiments in- 
volving point contacts at finite bias. Idealising the experimental situation and extending the notion of a quantum 
quench, we consider time evolution from an initial non-equilibrium state in a translationally invariant system. 
We show that electron interactions bring the system into a steady state at long times. Strikingly, this state is not a 
thermal one: its properties depend on the full functional form of the initial electron distribution, and not simply 
on the initial energy density. Further, we demonstrate that measurements of the tunneling density of states at 
long times can yield either an over-estimate or an under-estimate of the energy density, depending on details of 
the analysis, and discuss this finding in connection with an apparent energy loss observed experimentally. More 
specifically, we treat several separate cases: for filling factor v =1 we discuss relaxation due to finite-range or 
Coulomb interactions between electrons in the same channel, and for filling factor i/=2 we examine relaxation 
due to contact interactions between electrons in different channels. In both instances we calculate analytically 
the long-time asymptotics of the single-particle correlation function. These results are supported by an exact 
solution at arbitrary time for the problem of relaxation at ;/ = 2 from an initial state in which the two channels 
have electron distributions that are both thermal but with unequal temperatures, for which we also examine the 
tunneling density of states. 

PACS numbers; 71.10.Pm, 73.23.-b, 73.43.-f, 42.25.Hz 
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I. INTRODUCTION 

Recent experiments have stimulated renewed interest in the 
process of equilibration for an isolated quantum system. Some 
of the most prominent of these studies have involved time 
evolution in cold atomic gases, after a sudden change in sys- 
tem parameters.'-^ Other measurements, remarkably, are able 
to probe coherent quantum dynamics in conventional con- 
densed matter systems far from equilibrium. ^^ '^ Theoretical 
work has included studies of general features of dynamics 
after a quantum quench,"' and of the approach to thermal 
equilibrium," as well as detailed investigations of individual 
model systems»i2ii2 

A notable result from experiments on cold atomic gases 
is the absence of relaxation in some cases, despite there be- 
ing strong interactions between particlesji This has been at- 
tributed to the fact that the system concerned is, to a good ap- 
proximation, described by the Hamiltonian for an integrable 
model. 

In the condensed matter setting, a series of recent inves- 
tigations have examined aspects of the physics of quantum 
Hall edge states far from equilibrium. These include: the 
discovery of interaction effects in electronic Mach-Zehnder 
interferometers;^ the detection of thermal transport by quan- 
tum Hall edge states;"* and the measurement^ of the energy 
distribution of electrons in an edge state driven out of equi- 
librium by passage through a quantum point contact (QPC) at 
finite bias voltage, and the observation^ and control^ of the 
relaxation of this distribution. 

The theoretical work we describe in this paper is motivated 
in part by the last of these sets of experiments. In outline, 
the experiments^^ involve bringing together, at a QPC, in- 
teger quantum Hall edge states that originate from contacts 
at different potentials. Tunneling between edges within the 
QPC generates a non-equilibrium electron energy distribution 



downstream from the contact. This distribution is measured^ 
using a quantum dot as a spectrometer, as a function of propa- 
gation distance.^ It relaxes to a stationary form, which is close 
to thermal but with an effective temperature higher than that 
of the sample as a whole. The fact that this distribution propa- 
gates further with little change suggests that electrons in edge 
states are weakly coupled to other degrees of freedom. 

These observations raise some important theoretical ques- 
tions. A central one stems from the fact that the usual, chiral 
Luttinger liquid model for quantum Hall edge states''* is in- 
tegrable, being simply quadratic in the edge magnetoplasmon 
coordinates. Within this description, the number of excitation 
quanta in each collective mode is conserved under time evo- 
lution. In this context, it is natural to ask how relaxation takes 
place in edge states, and what steady state is reached at long 
times. 

The main difficulty in addressing these issues arises be- 
cause interactions are simplest to handle using bosonization, 
while a QPC is most naturally represented using fermionic 
variables. The two existing theoretical discussions'^-'^ of 
edge state relaxation aim at a faithful, but necessarily approx- 
imate, treatment of the experimental aiTangement, including 
the QPC. Our objective in this paper is to explore a comple- 
mentary direction. We study interaction effects on time evolu- 
tion in a translationally invariant system, omitting the QPC al- 
together from our model Hamiltonian. In its place, we choose 
the initial state of the system to have an electron distribution in 
momentum like that generated when non-interacting fermions 
pass through a tunneling contact at finite bias voltage. Our 
problem therefore has the features of a quantum quench,'*' ex- 
cept that the initial state is not chosen to be the ground state 
of a simple Hamiltonian. In our treatment we expect evolu- 
tion as a function of time to correspond qualitatively to evolu- 
tion in the experimental system as a function of distance from 
the QPC, with time and distance in the two pictures related 
by a characteristic edge state velocity. While this correspon- 
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dence is not precise, our exact treatment of interactions brings 
advantages. In particular, we are able to show that the one- 
electron correlation function reaches a stationary form in the 
long-time limit that is not thermal. It would probably be hard 
to know whether this conclusion was reliable if it were the re- 
sult of an approximate calculation. Our calculations draw on 
theoretical ideas developed in recent discussions'^ ''^ of non- 
equilibrium effects in Mach-Zehnder interferometers. 

The remainder of the paper is organised as follows. In Sec- 
tion we give an overview of the results and discuss the re- 
lation of our approach to experiments and to previous theory. 
In Section |III] we introduce the model and notation we use. 
In Section |IV] we show how plasmon dispersion arising from 
finite range or Coulomb interactions in an edge at v — 1 gen- 
erates a stationary electron distribution at long times, for an 
initial state that has no many-particle correlations but it is oth- 
erwise general. We then specialise this treatment to an initial 
state in which the electron distribution in momentum consists 
of a double step. In Section[V]we treat a system atv = 2 with 
contact interactions between electrons in different channels. 
We obtain an exact solution for the time evolution with edges 
initially at different temperatures, and study the long-time dis- 
tribution reached for a general initial state. Finally, in Section 
IVIIwe use this solved example to discuss the tunneling den- 
sity of states and implications for experimental estimates of 
the energy density in the steady state. 



II. OVERVIEW OF RESULTS 

In the experiments- we are concerned with, edge channels 
of a quantum Hall system at filling factor i/ = 2 are driven out 
of equilibrium by applying bias voltage F to a QPC. Edges 
that originate from contacts with different potentials meet at 
the QPC where, depending on its width, either the inner or the 
outer channels of the two edges are coupled by quantum tun- 
neling. As a result, immediately after the QPC, the electrons 
in these channels have a non-equilibrium energy distribution. 
For example, in a non-interacting system at zero temperature, 
this distribution has two steps with energy separation eV, as 
shown in Fig. [1] The energy distribution of electrons is ex- 
pected to evolve as a function of distance downstream from 
the QPC, because of inelastic scattering induced by interac- 
tions. This distribution is measured in the outer channel at 
successive distances from the QPC, using a quantum dot with 
a single active electronic level. It is found that the double-step 
distribution generated by the QPC equilibrates after a prop- 
agation length of order 10 /xm. The observed equilibrium 
distribution is similar to a Fermi function, with an effective 
temperature close to that set by the average electron energy 
density in the channels just after the QPC. 

Two different approaches have been used in previous work 
to model these experiments. Onoi^ is formulated in terms of 
electron operators, using a Boltzmann-like equation to study 
the evolution of the electron distribution with distance as a 
consequence of electron-electron interactions. This method 
has the advantage that the initial electron distribution can be 
calculated simply by using a single-particle treatment of the 



QPC. Its limitation stems from the fact that the integrability 
of the usual, chiral Luttinger liquid model for edge states is 
presumably not respected. The other approach'^ is formulated 
in terms of edge magnetoplasmons. This allows a detailed 
treatment of interactions, at the cost of a phenomenological 
treatment of the QPC, in which the plasmon distribution just 
after the QPC is estimated from shot noise in the charge. 

We evade some of these difficulties by treating relaxation 
of an initial distribution in a translation-invariant edge. The 
relaxation is caused by scattering of electrons between differ- 
ent single-particle momentum eigenstates. By contrast, the 
collective plasmon modes of the system do not relax, because 
within standard approximations the bosonized Hamiltonian is 
quadratic. Instead, using the language of plasmons, relaxation 
occurs because in the interacting system different modes have 
different group velocities. To expand on this viewpoint, we 
use it to estimate the relaxation time for an initial state in 
which electrons occupy single-particle levels independently, 
with a specified energy distribution. This state can be viewed 
in terms of particle and hole excitations above a Fermi sea. It 
is characterised by the mean separation Ljnit between these 
excitations. For an electron energy distribution having two 
steps at energy separation eV, with edge velocity v, we have 
iinit ^ hv/eV. Because of plasmon dispersion, the width 
I of particle and hole wavepackets broadens with time t: the 
relaxation time t is the time scale at which their width is com- 
parable to their separation. 

The simplest application of these ideas is provided by a sys- 
tem at filling factor v — 2 with contact interactions between 
electrons in different channels. Plasmon dispersion in this sit- 
uation is characterised by velocities u_ and w+ for slow modes 
and fast modes, as reviewed in Section |V] Then the width of 
the wavepacket grows as ^ ^ {v+ — V-)t and the relaxation 
time is 

V+ + V- h 

T= -■ (1) 

v+ — V- eV 

Note that this timescale and the corresponding lengthscale 
(similar to one identified^' in Ref. vary inversely with 
bias voltage V. It would be interesting to search for such a 
dependence in experiment: this might be done from the data 
in Fig. 2 of the Supplementary Material to Ref. 6, by extract- 
ing the decay length of Toxc and analysing its dependence on 
SVu. 

As a second example, consider a single edge channel at fill- 
ing factor 1/ = 1 with finite range interactions. The plasmon 
frequency Uq varies with wavevector q, having the small-q ex- 
pansion 

0Jq = vq-iv/b)ibqf/3..., (2) 

where b > is a length that characterises the range and 
strength of interactions. In place of v+ — w_ we take the 
difference between the group velocity dqUJq at q ~ and 
at q = which is v{b/e'f. Setting £ = vt{b/lY yields 
t = b{vt/bY^^ and from t = L\^\t we obtain 
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Thus at V = 1 with finite range interactions, we find a 
timescale and lengthscale that vary with bias vohage as V^^, 
in striking contrast to behaviour at = 2. Again, it would be 
of interest to test this in experiment. 

For unscreened Coulomb interactions the plasmon disper- 
sion is 



= q['v + u\n{l/bq)] 



(4) 



and the width of the electron wave-packet grows with time as 
£ = ut, giving 



hu 

evV 



(5) 



Similar scaling relations have been discussed in Ref.Ty'in the 
context of the decoherence in electronic Mach-Zehnder inter- 
ferometers. 

We show in the remainder of this paper that at times much 
longer than t, the equal time, one-electron correlation func- 
tion reaches a stationary form. For v = 1 and an initial state 
in which single particle levels are occupied independently, we 
determine this form exactly as a functional of the initial dis- 
tribution. Its short-distance behaviour is fixed by the energy 
density in the nonequilibrium state, and characterised by an 
effective temperature. Its behaviour at long distances, how- 
ever, depends on the full functional form of the initial distribu- 
tion. We present detailed results for the case of a double-step 
initial distribution, as produced at a QPC. For this example 
we also consider the electron momentum distribution at long 
times, which shows power-law behavior at small momenta, 
with an exponent which depends on tunneling probability of 
the QPC. At large momenta the distribution is close to thermal 
with an effective temperature given by the excess energy. 

The electron distribution is probed in the experiments of 
Refs. Jj-tS, by measurement of the tunneling current through 
a quantum dot that has a single active energy level weakly 
coupled to the edge state. Without interactions, this current 
is proportional to the occupation probability of edge-state or- 
bitals at the energy level of the dot, and so yields a straight- 
forward determination of the electron distribution in energy. 
With interactions, however, the situation is more complicated, 
since matrix elements for tunneling between the edge and the 
dot depend on the many-body state of the quantum Hall edge. 
We examine the consequences of interaction effects for the 
interpretation of tunneling current measurements in Section 
IVII focussing in particular on estimates of the energy den- 
sity. We find that the energy density can be over-estimated or 
under-estimated, depending on how such measurements are 
analysed. 



III. MODEL AND OBSERVABLES 

In the following we study edge states in the integer quan- 
tum Hall regime at filling factors v = \ and v = 2. The 
microscopic Hamiltonian is a sum 



of contributions from kinetic energy Hkin and interaction en- 
ergy Hint- For Uc chiral copropagating electron channels, 
with channel index and Fermi-velocities u,,, we have 



rig „ 
ihJ^Vr^ J 



L/2 

-in y v„ I tJ)^ {x)dxiprjix)dx . (7) 



Here we impose periodic boundary conditions on a system of 
length L so that allowed wavevectors are fc = 2TTn/L, with in- 
teger n. Fermionic creation and annihilation operators c^^^ and 
Ckrj for an electron with wavevector k on the edge 77 obey stan- 
dard anticommutation relations {ckr),Cpjj,} = SkpSjfjff. The 

field operator iprjix), which annihilates an electron at position 
X on the channel rj, is 



k— — oo 



ikx 



(8) 



We consider translationally-invariant interactions described 
by the potential {x — x'), so that 

I rL/2 ,L/2 
Hint = - / dxdx' Urin'ix - x')pr,{x)pjj>{x'). 

^ J -L/2 J -L/2 

(9) 

We discuss two different cases; a single edge channel at 
V = 1, with finite range interactions; and a two-channel sys- 
tem at = 2, with contact interactions between electrons in 
different channels. Contact interactions between electrons in 
the same channel can be absorbed into a renormalization of 
Fermi velocities. 

Our first objective is to obtain the equal time single-particle 
correlation function. Using a time argument to indicate oper- 
ators 'iljjj{x,t) = e*^*?/',,(x)e~*^* in their Heisenberg repre- 
sentation, we write this quantity at time t as 



G„{x-X',t) = {lP+{x,t)i^r,{x',t)), 



(10) 



where the average is calculated with respect to the initial state, 
yet to be specified. The electron momentum distribution in a 
translationally invariant system is then 



L/2 

n^{k,t) = I dx e'^'^Grfix.t). 

-L/2 



(11) 



In thermal equilibrium at inverse temperature /3 the single- 
particle correlation function has at = 1 the form 



G{x, t) = i {2/3;iusinh[7r(a; - vt + ia)/l3hv\y 



(12) 



with a lengthscale fihv, timescale f3h and a short-distance cut- 
off a. 

The experiments of Refs. IH-ItI measure the current through 
a quantum dot (QD) weakly coupled to one edge channel. As- 
suming sequential incoherent tunneling through a QD with a 
single energy level E, perturbation theory in the tunneling am- 
plitude tr) gives for this current the expressions^ 



H ^ H, 



kin 



Hi. 



(6) 



e\tD? 



(13) 
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where the tunnehng density of states is at long times 

/oo 



IV. EQUILIBRATION Xt v=l 

To compute the time evolution of operators under H we 
use bosonization^^ employing the same notation as in section 
IVC of Ref.[T^ Consider a system at = 1. The Hamiltonian 
in bosonic form is 

H = J2f^ib+b<i (15) 

q>0 

where 6+ is an operator creating a plasmon with wavevector 
q and energy huq, and we omit the channel label rj. Let the 
Fourier transform of the interaction potential be 

/oo 
dx e-^«^C/(x) . (16) 
-oo 

Then the boson dispersion relation is 

u!q = vq[l + u{q)/2Trhv] . 

In this section we express the correlation function, Eq. ( fTOb . 
in terms of the electron density operator p{x) — ip^ {x)tp{x) 
written in the Schrodinger representation. The result, Eq. ( fT9l ), 
involves an average of p{y) with a kernel, denoted below by 
K{x, t; y). Quantum fluctuations of this average are Gaussian 
at long times, because the range of the kernel becomes much 
larger than the correlation length in the initial state. For this 
reason, the correlation function in the long time limit can be 
expressed as a functional of the second moment of the density 
in the initial state. More formally, one can write a cumulant 
expansion for the correlation function, in which the second 
order cumulant is time independent and higher order contri- 
butions decay to zero at long times if the plasmon dispersion 
Ljq is not exactly linear in q. 

A. Correlations at = 1 in terms of initial density operator 
Introduction of the bosonic field 

q>0 

allows the electron correlation function to be written as 

G{x,t) ^ (27ra)-i(e''^(^'*'e-*'^(°^*)) . (17) 

In order to evaluate this average, we express the boson opera- 
tors bq and 6+ in terms of p{x). From the standard expression 

P^y) = ^ J2i^7r/qLy/^[iqbqe^'^y - iqb+e-^^y], 

q>0 



we have for fc > 

/■i/2 r /o \l/2 

Thus the difference in bosonic fields 

rL/2 

(j>{x,t)-<l){Q,t)^ K{x,t-y)p{v)dy, 

J-L/2 

is given by a spatial average of the electron density operator 
weighted with the kernel 

K{x,t;y) = -iXl^(e-*9^ - l)e*'^^*+*'3^ (18) 

This yields an expression for the correlation function in terms 
of fermionic operators 

G{x,t) = {2Tra)-^e"^'=''''>/^{e' I ^^'='*'y'>f''^y'>) , (19) 
where 

I{x, a) EE [<P{x, t), 0(0, t)] = ln[(a + ix)/{a - ix)]. (20) 

1. Long-time asymptotics of the kernel 

Considered as a function of y, the kernel K {x, t; y) at time 
t = is a step of width \x\ with height —2tt. In the non- 
interacting system for t > {), this step propagates at velocity 
V without change in shape. With interactions, the relaxation 
time of r introduced in Section II is the timescale for the step 
to change shape. We are interested in behavior of K at long 
times in presence of interactions. Consider first finite-range 
interactions with the plasmon dispersion at small wavevector 
given by Eq. (|2]l. We separate out the average motion by writ- 
ing y — x/2 — vt + X. This gives for the phase appearing in 
Eq. (HI 

q{y - x/2) + LOqt = qX- {vt/b){bqf/i + ... (21) 

at small q. Defining a lengthscale Is = b{vt/bY/'^ as in Sec- 
tion |II] the kernel at long times has the form 

K{x, t, y) ^^2 r ^ sin ^ cos[Q{X/ls) + Q'/S] . 

(22) 

Explicit evaluation for \x\ ^ Is yields 

Kix, t, y) = ~2ti(x /ls)ki[X/ls] • (23) 

Thus the range of the kernel grows with time as t^/^ and its 
amplitude decays as t^^/^ . 

We can obtain asymptotics of the kernel for a system with 
Coulomb interactions in a similar way. In this case the plas- 
mon dispersion is given by Eq. (HJi. Writing y — x/2 — vt — 
uthi{ut/b) + ^, the phase appearing in Eq. ( fTsT i is 

q{y - x/2) + ujqt = gA + qutln.{l/qut) . (24) 
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Defining Ic — ut, we have 



Kix, t, y) = -2 r ^ sin ^ cosiQX/lc - Q In |Q|)] , 

(25) 

Hence, with Coulomb interactions at long time, the range of 
the kernel grows as t and its amplitude decays ast^^. 



2. Cumulant expansion 

To evaluate the average in Eq. ( fT9] l we write it in terms of 
cumulants, as 



exp 



n i-L/2 



L/2 



dy dzi . . . dzn-i 



X Fn{x,t,y;zi . . .z„_i)C„ (zi . . .z„_i) 



(26) 



where n-th order cumulant of the electron density p{y) is 
C„ (zi, ...z„_i) = {p (y) p {y + zi) . . . p {y + Zn-i)) (27) 
and the rt-th order convolution of kernels is 
1 



Fn{x,t; Zi . . . Zn-l) 



2tt 



dy K{x,t;y) 



-L/2 



X K{x, t;y + zi) . . .K {x,t;y + z„_i) 
Consider the second order contribution. This is 

rL/2 



1 

F2{x,t;z) = — dy K{x, t; y)K{x, t 

J-L/2 



y + z) 



L/2 



dy e 



i(q+k)y ikz i(u}q+u}^)t 



(28) 



-L/2 



The result, independent of both time and interactions, is 

F2{x,z)^Ti{\x + z\ + \x-z\-2\z\). (29) 

It is a triangle as a function of z, which is non-zero in the 
region (— x) and has amplitude ^2(2^, 0, 0) = 2ti \x\. 
At third order we have 

/ ,3 dq f°° dp 1 
F^{x,t-z^,Z2) = -{-if ^ 

J -00 q J -00 p p + q 

X (e-"^'= - 1) (e-'P'' - 1) (e^(P+9)^ - 1) 

X g«('^p+'^<!-'^p+<!)*gip2;ig-i(p+(/)z2 

Time enters the integrand of this expression through the factor 
gj(wp+w,-wp+,)t^ Without interactions ujq + ujp - ujp+q = 0, 
but in an interacting system this factor oscillates rapidly at 
long times, as a function of p and q. For large t the inte- 
grals are therefore small and are dominated by low-frequency 



contributions. In this regime we can use the expansion 
+ — ijjp+q ex. q'^ + p^ — [p + q]'^ and remove the time 
dependence from the exponent by rescaling variables, with 
q = Q/ls^ P — P/ls, as in the discussion of the kernel, and 
X = Xls, Z1.2 = ^i,2^s ■ In this way we find F3 = i^^/'^x (a 
time-independent function of X, Zi, Z2). Because the func- 
tions Cn (zi . . . Zn-l) decay at large z we can use in Eq. ( l26b 
the function F3 evaluated at Z\ — Z2 — 0, obtaining a contri- 
bution that decays as t^^l'^ . An equivalent result holds for all 
higher order cumulants. Similar arguments also apply in the 
case of Coulomb interactions, but with third and higher order 
terms decaying in time as 

The significance of these results is that contributions from 
all cumulants beyond the second one decay with time and 
may be neglected in the long-time limit. Cancellation of 
the high order cumulants is known for thermal states as the 
Dzyaloshinskii-Larkin theorem, '^'^"^ but is not in general valid 
for nonequilibrium states. We have thus shown that a similar 
result holds in the long-time limit. 

In summary, we obtain for the correlation function in the 
long time limit the expression 



Q[x) = lim G(x,t) = (27ra)^ie^^('')x 



exp 



— TT / dz F2{x,z){p{x)p{x ^ z)) 



(30) 



This is one of the main results of this paper. In order to use it, 
we need to calculate {p{y)p{y + z)) for the initial state, which 
we do below. 



B. Long-time correlations 

As a check and because we shall need to refer to the results, 
we start by considering averages in the ground state, denoted 
by (. . .)o. We require (/5(0)/5(z))o and find 



(/5(0)p(z))o = 



1 



1 



(27r)2 (a + ^z)2■ 
Integrating this with the function F2 {x, z) [Eq. ( [29l )1 we have 

dz F2{x, z){p{x)p{x + z))o 



= - — [Ina^ -ln{x^ + a^)]. (31) 

Substituting Eqns. ( ISTT i and ( l20t into Eq. ( l30t we obtain for 
the ground state correlation function 



G{x,0) = (2Tra)-^exp{^[\na^ - \n{x^ + a^)]} 



X [{a + ix)/{a — ix)] 



1/2 _ 



1 



2Tr X + ia ' 



in agreement with Eq. ( fT2] i. evaluated for P 00. 

Now consider an initial state in which single-particle mo- 
mentum eigenstates are occupied independently with prob- 
ability n(fc) as a function of wavevector k. We separate 
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this probability into its value no(fc) for a filled Fermi sea 
with Fermi wavevector zero, and a deviation 6n{k), writing 
n{k) = no(fc) + 6n{k). Without loss of generality, we take 



Sn{k) dk = 0. 

In this state we have 

fclfc2fc3 

= (p(0)p(z))o + C(^'(z) + cf (z) 
where we have introduced 



(32) 



2)Z 



and 



1 



dqdn{q) sinqz . 



1 

2^ 



dq 5n{q)e-"i' 



So in the long time limit we have 

: 1 



g{x) 



2tt X + ia 



cxp[-/ (a;) 



(33) 



(34) 



(35) 



(36) 



with 



fix)^7r dzF2{x,z)[ci'\z) + ci^\z)]. (37) 



This represents a second key result of our work, and we next 
analyse the behaviour of f{x). 



1. Short-distance asymptotics 

Consider f{x) for x much smaller than the quasiparticle 
spacing in the initial state. Since F{x, z) is non-zero only 
for \z\ < \x\ we require the small z behavior of C2^\z) and 

(2) 

C^ '{z). Both can be expressed in terms of the energy density 
in the initial state, characterised by an effective temperature 
T*. This is defined following the relation 

^ j 6n{q)qdq=^{k^Tf (38) 

between energy density and temperature T for a thermal 
state. We have ^(z) = {ksT* / hv p)^ / 12 + ©(z^) and 
^(z) = -[TTzCf In this way we find 



fix) 



7r2 //crT*" ^ 



6 V hvp 



x^ + Oix^) . 



(39) 



In summary, the short-distance behaviour of the steady state is 
identical to that in thermal equilibrium at a temperature fixed 
by the initial energy density. 



2. Long-distance asymptotics 

At large x, the function F2 (x, z) has slow dependence on z, 
but c'^^ (z) and Cj^'' (z) fall to zero for |z| much greater than 
the quasiparticle spacing in the initial state. We can therefore 
make the replacement F2{x, z) = 2n\x\ in Eq. dJTl l. obtaining 

/OO 
dq [sgn{q)5n{q)~\5n{q)\'] + 0{\x\^) . 
-00 

(40) 

The long-distance decay of the correlation function is there- 
fore not in general determined solely by the energy density, 
but instead depends on the entire momentum distribution in 
the initial state. In the following subsection we examine the 
consequences of these general results for a specific choice of 
initial distribution. 



3. Equilibration from double-step initial momentum distribution 

We discuss in detail the long-time correlations for an initial 
state that has a double-step momentum distribution as shown 
in the Fig. [T] which serves as a model for the electron distri- 
bution generated just after a QPC in tunneling experiments.*^ 
For a QPC with a tunneling probability p at bias voltage V, 
steps in 5n{q) are located aX q = —Q- and q = Q+, with 
separation + = Q = eV/hvp, so that 



6n{q) 



-{l-p) -Q-<q<0 
p <q <Q+ 



and Sn{q) = if q lies outside this range. Eq. (|32] | implies 
(1 — p)Q- = pQ+- With these definitions we have an effec- 
tive temperature 



T* = V3 {hvFhks) QVp{1 - P) ■ 



(41) 



It is interesting to examine for this initial distribution how 
the steady-state correlation function differs at long distances 
from a thermal one with the same effective temperature. We 
have 



dq\5n{q)\' =Qp{l-p), 



' n(q) 
















P 

q 




Q 





FIG. 1: Double step electron momentum distribution. The width Q 
of the step represents a voltage hvQ applied to a QPC, which has 
tunneling probability p. The quantities Q+ and Q_ are determined 
following Eq. ([32}. 
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and 

/oo 
dq sgn{q) Sn{q) ^ 2Qp{l - p) 
-oo 

so that for (3|a;| > 1 

f{x)=np{l-p)Q\x\. (42) 
For comparison, in a thermal state one has from Eq. (fT2l l 



\x\. 



(43) 



A simple replacement of T with T*, expressed in terms of 
p and Q using Eq. (HTt . would lead one to expect in place of 
Eq. (|42li the result (also fovQ\x\ > 1) 



(44) 



Clearly, (except for the special value p(l — p) = Stt^'') 
these are distinct forms. The difference is most pronounced 
if p{l — p) ^ \, when the correlation function in the steady 
state decays much more rapidly with |x| than would be the 
case in a thermal state with the same energy density. 

For this case of a double-step initial distribution, the steady- 
state correlation function for general x can be obtained explic- 
itly, in the form 



g{x) 



2tt X + ia 
where we have defined the function 



e^^[-p{l-p)R{Qx)l (45) 



1 



II 



- 2(1 + 7£; + Inz - cosz - Ciz - zSiz), (46) 

Si and Ci are the sine and cosine integral functions, and ~ 
0.577216 is the Euler constant. 



4. Resulting long-time momentum distribution 

It is also interesting to discuss the steady-state momentum 
distribution within this framework, especially in the regime 
p{l — p) ^ \ where deviations from thermal behaviour are 
largest. As a first step we use Eq. ( l45T l to derive a form for 
the real-space correlation function at large \x\ that includes 
sub-leading terms omitted from Eq. (|37] ). obtaining 



g{x) 



-7rp(l-p)|Qa:| 



27r 



X + ia 



(47) 



From this, using Eq. (fTTT i. we find that the stationary momen- 
tum distribution for |fc| ^ Q is 



n{k, t) 



1 1 

2~ n 



ga(i+7i=)r(Q,)lni[(-i/c/Q + 7ra/2)-"] , 

(48) 



where a = 2p{l~p). A striking point is that two scales appear 
in Eq. (|48] |: one is Q, while the second one is ■KaQ/2. These 
scales are widely separated when a is small, and so there are 
two distinct small k regimes of behavior for n{k,t) at long 
times, with 



n{k,t) = i - C(7ra/2)~(i+")afc/Q 



(49) 



for k <^ TTaQ/2 and 



1 



n{k, 2 ~ ^(W2)[l - aQ/k]{Qlk)" (50) 

for 7rag/2 k <^ Q, where C = ie"(i+T^)r(a). 

As a supplement to this discussion of asymptotic forms, the 
full momentum distribution at long times can be obtained nu- 
merically from Eqns. (l45b and (fTTl l. We present its behaviour 
in Fig. |2] Since the quantity measured experimentally is a 
derivative, we plot the difference in derivatives between the 
steady state distribution and the one for a thermal state with 
the same energy density. As anticipated in our earUer discus- 
sion, these differences are largest when p{l — p) is small. 



V. EQUILIBRATION AT = 2 

In this section we study relaxation of a nonequilibrium dis- 
tribution in a system of two channels coupled by contact inter- 
actions. Whereas contact interactions at z^ = 1 simply renor- 
malise the Fermi velocity, without providing a mechanism for 
relaxation, at z^ = 2 they give rise to two linearly dispersing 
collective modes with distinct velocities, and hence do pro- 
duce relaxation. We therefore focus on this as the simplest 
interesting case. We discuss briefly the effects of finite range 
interactions at = 2 in Section [VC3l 



0.05 - 



-0.05 - 



O 



-0.15 



-0.25 




FIG. 2: Difference between derivatives of the electron momentum 
distribution functions in an steady state described by Eq. ( 145 1 and 
in a thermal state with the same energy density. Thick solid line: 
p = 0.5; dotted line: p — 0.3; short-dashed line: p — 0.25; long- 
dashed line: p = 0.2; thin solid line: p — 0.15; dot-dashed line: 
p = 0.1. 
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The Hamiltonian of the system is given by Eq. (|6]l, with 
ric — 2 and an interchannel interaction potential 

Ui2{x - x') = U2i{x - x') ^ {g/2)6{x - x') . 

Intrachannel interactions are absorbed into the Fermi veloci- 
ties Wi.2- 

We start this section with a general discussion of time evo- 
lution for this Hamiltonian. Then we set out the solution 
to a special case, calculating the electron correlation func- 
tion at arbitrary time for a particular choice of initial state, in 
which each channel is separately in thermal equilibrium, but 
at different temperatures. Finally, we extend the approach de- 
scribed in Section |IV] to calculate the long-times asymptotics 
of the correlation function at long and short distances for an 
arbitrary initial single-particle distribution. 



using the transformation S and its inverse. In this way we 
obtain 



ak{t) = Uk{t)ak +Sk{t)bk 
and bk{t) = Sk{t)ak + Ok{t)bk. 

with coefficients 



Ufc(t) = 

_^ — iv-\-kt ^ — iv—kt\ 



u-kt -2 



cos^ 



(56) 



sin^ e, (57) 



and Sk{t) — [e — e '""""'j cost/smf 

Equations (|56] | and (ISTT i define a unitary transformation of the 
boson operators which generates their time evolution indepen- 
dently for each wavevector. 

In bosonized form, the electron correlation function is 



A. Diagonalization of bosonic Hamiltonian 

The bosonized form of the Hamiltonian is 

H — hvi ka^hk + hv2 kb^bk 



k>0 



k>0 



\^T.^i^^^t + ^4bk), (51) 



2 27r 



fc>o 



where Ofc and bk are the boson operators on the first and the 
second channel. It can be diagonalized by the rotation 



{ak.Pkf = Siak,bk)'^ , 



(52) 



with 



„ _ , cost/ smf, ■ 
^" ' -sine cosP I ■ ^^^^ 



Setting 7 = 5/ (27r/i) and tan 29 — 7/(^1 — f 2), this brings 
the Hamiltonian to the diagonal form 

H = hv+Y, ka+ttk + hv- kp+Pk ■ (54) 

fc>0 fc>0 

The collective mode velocities v± are 

v+ —vi cos^ 6 + V2 sin^ ^ + ^7 sin 29, 
u_ =ui sin^ 9 + V2 cos^ - \l sin 29. 

The time dependence of the operators dk and Pk in the 
Heisenberg representation is 

= e-'"+'^-*afe , Pk{t) = e-™-'=*/3fc • (55) 

Since we wish to discuss time evolution from an initial 
state specified by occupation numbers for independent, non- 
interacting fermions, we need to re-express Eq. (l55l l in terms 
of the operators a and b in their Schrodinger representation. 



G,,{x,t) ^ {e"'i;+{x)^^^{Q)e-^"') 

= (27ra)-i(e'''^'''=''*^e-*"^''(°'*^) . (58) 

For an initial state without correlations between the channels, 
this factorises into the product 



G^{x,t) = {2^a)-^G^:^\x,t)G^^\x,t), (59) 



where 



and [from Eq. ( l56] ll 

4>^^\x,t) ^-Y,{2i,lqL)'/\n,{t)a,e^'>- +h.c.)e-'"^l\ 

q>a 

(60) 

The corresponding expression for 0^^^ (x, t) is obtained by 
making the substitutions Uq{t) — !> 5q{t) and — > bq in Eq. 
(I6OI 1. Similarly, we obtain an expression for bosonic fields in 
the second channel by making the substitutions Sq Vq and 

Uq — ?> Sq in the equations for 4>[^l{x, t). 



B. Equilibration of channels with unequal initial temperatures 

We next study relaxation from an initial state in which the 
particles in channels of the system without interactions have 
thermal distributions with different temperatures, T^,. A simi- 
lar problem was studied in a context of a quantum quench for 
a non-chiral Luttinger liquid, in Ref . [25l 

The simplifying feature of a thermal initial state is that there 
is uncorrected occupation, both of the fermion orbitals on 
which Ck,-i act, and of the boson orbitals on which and bk 
act. We can therefore specify this state by the averages 

(a+afe) = (e'"'i'=^i-l)-\ = (e'"'^'''^-l)"\ (61) 

where = l/fc^T^. Correlators of the fields 4>]^\{x,t) are 
those for a non-interacting edge, with the appropriate choice 
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for velocity. We use the notation Wr, = —2ij3nhvri/ L and 
= n/Pj^hvrj. We also define r = ^ sin^ 29, s = 1 — r, and 



vd = v+ 



7 



sin 26* ' 



(62) 



In this way we obtain 



1 



1 



wf sinh*[^i(x + ia)] 
sinh[^i(w£)i + ia)] sinh[^i(— w^t + ia)] 



sinh[^i(a; + vjjt + ia)] sinh[^i(x — vut + ia)] 



In a similar way we get for G^^^ {x, t) 
G^^\x,t)^{2na/Lr ^ 



r/2 



(63) 



W2 sinh''[^2(a; + ia)] 



sinh[^2(a; + vot + ia)] sinh[^2(-T — vot + ia)] ^^"^ 



(64) 



sinh[^2('UD^ + ia)] sinh[^2(— + ia)] 
As a check, note that the initial electron correlation function 

(G^(a:,0)) = (z/2^)C^sinh-i[e^(2; + m)] 

derived in this way from Eq. ( l59b coincides with the thermal 
one at inverse temperature Next consider the long-time 
asymptotics of Grj{x,t). The function 1/ sinh[^j,(a; + ia)] is 
peaked near x = and falls off exponentially at large x on a 
scale 1 / (or as a power-law at zero temperature). For a fixed 
X at long times, the factor [. . .]''/^ approaches one. The length- 
scale defined by the product sinh[^i(a:: + ia)]~^ sinh[^2(a: + 
ia)]~'^ is 



(sCi+r6)-' 



PlVi I32V2 



The time taken to approach the asymptotic form is t l/vjj. 
For times much longer than r we have 



C. Correlations at i'—2 in terms of initial density operator 

As a second example of relaxation in the two-channel sys- 
tem, we consider an initial state with a double step distribu- 
tion of electron momentum in one channel. The correlation 
function decouples into a product [Eq. ( |59] l1 of factors that 
can each be expressed using the approach of Section ITV Al in 
terms of the initial electron density operator. For example, 

G*,^'(a;,i) = e^^'^'^^^'^^''^'^'^^' '*^^e^ ■I' '*'^^^''^^'^^) , (66) 
where the kernel is given by 

K{x,x',t;y) = 5](27r/gL)4(0(e"*«^-e-'''^')e^'^-'^'" 

(67) 

and has the form 

K{x, x', t- y) = 7r( [sgn(w+t - x + y) 

— sgn(w+t — x' + y)] cos^ 9 
+ [sgn(w_t — X + y) — sgn(u_t — x' + y)] sin^ 9) . 

Attimesi > this kernel separates into two steps: 

one occupying the interval (a:' — v+t, x — v+t), and the other 
in {x' — V-t,x — v-t). The separation between these intervals 
grows at speed vot- An energy scale eV in the initial state 
sets a characteristic length hv /eV: the timescale r introduced 
in Sec. II is the time taken for the separation between steps to 
reach this characteristic length. Introducing an operator 



N{x) 



p{z)dz 



that counts the number of particles in an interval of length x, 
in the long time limit there is the factorization 



pi / K{x,0,t;y)p{y)dy\ 



-2i-!vN{x) E 



2h {(3ivi)' sinh'lS^iix + ia)] 



{P2V2Y sinh''[^2(a^ + ia)] 

We see that at long times the system reaches a steady state. 
In this state the electron correlation function is strikingly dif- 
ferent from a thermal one, being given by a product of two 
thermal correlation functions for systems at different temper- 
atures, each raised to powers that depend on the strength of the 
interchannel interactions. Comparison of electron momentum 
distribution obtained from Eq. (l65T l with a thermal one is pre- 
sented in Fig. [3] 

Clearly, an attraction of studying evolution from this spe- 
cial class of initial states is that one obtains explicit results at 
not only at long times but at all intermediate times. In Fig.|4] 
we exploit this to show the time evolution of the momentum 
distribution in a channel initially at zero temperature, coupled 
to a second channel with non-zero initial temperature. 



0.02 




-0.04 i ' ' ' • ' • ' • ' ' ' 1 

-6 -4 -2 2 4 6 

hvpkfi 



FIG. 3: Difference between derivatives of the electron momentum 
distribution functions in an steady state described by Eq. l i65t with 
T2 — and in a thermal state at temperature T = ^/sTi . Thick 
solid hne: 6 — 7r/4; dashed line: 9 = 7r/6; dot-dashed line: 9 = 
7r/8; dotted line: 9 = tv/12. 
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Here we see a contrast between our earlier results atu = 1 
in systems with finite range or Coulomb interactions, for 
which Uq has dispersion, and these results at = 2 with 
contact interactions between electrons in different channels, 
for which the bosonic modes have linear dispersion relations. 
While in the former case the steady-state correlation function 
can be expressed solely in terms of the two-point density cor- 
relation function of the initial state [see Eq. (I30]l1. this is not 
so in the latter case. Nevertheless, simplifications arise in the 
small and large x asymptotics of the correlation function, as 
we now show. 



operator is a well-studied problem. The asymptotics can be 
obtained using the theory of Toeplitz determinants: see for 
example Ref. 13. 

In the case of a double-step distribution of width Q, which 
is created and measured in channel 1, with channel 2 initially 
at zero temperature, we obtain for Q\x\ ^ 1 in the leading 
order 



Gi(x) 



1 



2tt X 



■ exp[ 



, \Qx\ 
■ 2n 



(69) 



where 



1. Short-distance asymptotics 

The short-distance form of the steady-state correlation 
function is fixed by the energy density in each channel in the 
initial state, characterised by effective temperatures T*. Set- 
ting S,* = ■nkBT*/hvn, we find 

Correlations at short distances in channel 1 therefore have an 
effective temperature 



s = --[ln(l - Ap{l-p) sin2[7rcos2 9]) 



*2 



(68) 



with the corresponding result for channel 2 obtained by ex- 
changing r and s. We see that the effective temperatures of 
the two channels are different, except in the strong-coupling 
limit where r = s = 1/2. 



2. Large-distance asymptotics 

Calculation of the large x asymptotics of the Green func- 
tion in the form given by exponent of the electron counting 





\ \ 




0.8 


\ 




0.6 






0.4 






0.2 




\ 


n 







FIG. 4: Time evolution of the electron momentum distribution for 
a system of two channels in the strong-couping limit (6* — 7r/4) 
with unequal initial temperatures Tj* = and T2 , obtained using 
Eqns. (59), (63), and (64). Dotted line: ^2VDt = 1; dashed line: 
^2VDt = 2; dot-dashed line: £,2VDt = 3. Solid line: thermal distri- 
bution at temperature T — T2 /\/2, shown for comparison. 



ln(l - 4p(l -p) sin2[7rsin2 6])]; 



see Refim^ for details. 



3. Finite range interactions at v — 2 

Finite range interactions, in contrast to contact interactions, 
generate dispersion so that the mode velocities v± and the 
mixing angle 6 become functions of wavevector. Dispersion 
in the two mode velocities provides a further mechanism for 
relaxation at = 2. One sees, however, from Eqns. (1) and 
(2) that this additional mechanism is less important at small V 
than the one due to contact interactions, since it yields a much 
longer relaxation time (scaling as V^^ in place ofV^^). Thus 
at small V there is window between these two time scales dur- 
ing which the results for contact interactions apply even in the 
presence of dispersion. At times beyond the longer scale, our 
results for = 1 carry over with simple modifications and we 
obtain for < 00 



Gl{x) = [ga{x)f 



9+sin'' ( 



2 cos^ e sin^ e 



(70) 



Here Qa{x) and Qb{x) are obtained from the expression for 
Q {x) in Eq. dSOl l by computing expectations values in the ini- 
tial stares of channel 1 and 2 respectively, while 6 denotes the 
limiting value of the mixing angle at small vectors. 



VI. TUNNELING DENSITY OF STATES 

As summarised in Section |lll the experiments of Refs. 111-0 
probe the electron distribution by means of tunneling mea- 
surements. With that as a motivation, we study in this sec- 
tion the tunneling density of states (TDOS). We consider in 
particular edge channels at = 2 after long-time evolution 
from an initial state characterised by two temperatures, as in- 
troduced in Section lVB] We focus on this case, rather than an 
initial state with a double-step distribution which would bet- 
ter represent the experimental situation, because the analytic 
results available for the Green function with these initial con- 
ditions facilitate accurate evaluation of the TDOS. It is im- 
portant to stress that, while for a non-interacting system the 
electron momentum distribution and the tunneling density of 



states are simply related, with interactions they become inde- 
pendent quantities. This is illustrated (in the space and time, 
rather than momentum and energy domains) in our results be- 
low by contrasting Eq. i65[ with Eq. ( l73b . 

For a single edge state in equilibrium at temperature T with 
only contact interactions, the TDOS defined in Eq. (fT4l l has 
the form 

i^iE) = {hv)-^nFiE/kBT), (71) 

and so the tunneling current is proportional to the Fermi dis- 
tribution npix) — (e^ + More generally, one can at- 
tempt to characterise the energy dependence of the TDOS by 
an effective temperature Toff- As there is no unique way to do 
so, we explore two alternative. First, since the energy density 
in thermal equilibrium is j^^ik-BT)^, one natural definition of 
Tcff is via 

[u{E) - u,ae{-E)]EdE , 
(72) 

where 9{E) is the Heaviside step function, and the normalisa- 
tion is j/cff = oo), the TDOS far from the Fermi energy. 
To evaluate v{E) we require [see Eq. (fT4l) 1 the correlator 

GJt)= \mi{i!+{x,t + T)'4>{x,t)). 

Consider the two-channel problem discussed in Section FVl 
In thermal equilibrium the correlator Gi (i) reads 



1 

^ sinh[^(-i;_t + ia)Y''^^ ^ ' 
which gives for the density of states 

v{E) = {hv''°''' '^y^npiE /keT). 

Thus the equilibrium TDOS in a two channel system with con- 
tact interactions between electrons in each channel is propor- 
tional to the Fermi distribution and independent of interaction 
strength. It is worth noting that this is not the case for a single 
channel with finite-range or Coulomb interactions. 

Now let us consider the tunneling density of states at long 
times, when the system starts out of equilibrium. In order to 
have explicit results that are easy to evaluate accurately for all 
energies, we consider an initial state with channels at unequal 
temperatures T*2 - For this case the Green function in channel 
"1" at long times reads 

Gi{t) = L-^w^'w^'' X 

1 



sinh[^2(-w+t + ia)Y/^ smh[^2{-V-t + m)]''/^ 
1 



sinh[^i(-u+t + m)]™'^'' " sinh[^i(-i;_i + ia)^''^'' ^ ' 

(73) 
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FIG. 5: Derivative of the normalized tunneling density of states as 
a function of energy in a two-channel system. Filled circles: in the 
steady state reached with initial temperatures Ti > and T2 = at 
g/hvi = 1 in the strong coupling limit vi — V2. Dashed line: in a 
thermal state at the temperature Ti /\/2 expected from equipartition. 
Solid line: in a thermal state with a temperature 0.64 x Tj* chosen to 
give the best fit to the filled circles. 

This gives for the integral Eq. (|72] | 

+ [vl + v^_){s + {viT;/v2T*fr}] . (74) 

The effective temperature T^s obtained from Eq. ( |74| ) is a 
function of two parameters: 9 and V2/V1. For vi — V2 and 
with the second channel at an initial temperature T2 = we 
obtain 

^(fcBTeff)^ = ^ikBT*ni + h/2v,n (75) 

From this we find Teff = {T^/V2)[l + (7/2t;i)2]i/2^ which 
is always higher than the equipartition value T* / a/2. 

An alternative way to define Tcff is to fit the normalized 
TDOS to a thermal Fermi function, as displayed in Fig. |5] 
Interestingly, with vi = V2 and for g/Hvi < 1, the best fit 
yields Tcs ^ 0.64 x T^, which is about 9% lower than the 
equipartition result. 

We believe these calculations of effective temperature make 
some useful points. First, they show that there is no unique 
definition of effective temperature because the steady state is 
not thermal. Second, they serve to demonstrate that it is possi- 
ble for effective temperature (at least, by our second definition 
and for the initial state we have treated) to be lower than ex- 
pected from equipartition. Further, the specific mechanism for 
a low effective temperature is clear: the derivative of the tun- 
neling density of states at long times has less weight in its peak 
and more weight in its flanks than in thermal equilibrium. 

In the experiments of Refs. 5 and 6 an apparent energy loss 
is reported from analysis of the TDOS far from the QPC, com- 
pared to behaviour close to the QPC. The experimental anal- 
ysis uses two methods to find an effective temperature, corre- 
sponding to the two definitions we have employed. Although 
no systematic discrepancies are reported between results from 
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the two methods, the first method (calculation of energy den- 
sity by integration) relies on capturing contributions from the 
flanks of the distribution, which in turn requires an accurate 
value for the baseline. It may therefore be less robust than the 
second method (fitting to a Fermi function). We find it inter- 
esting to note that the temperature reduction we obtain from 
this method is within the errorbars of the found in Ref. i6|. 
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